Casimir-Polder interaction between a polarizable particle and a plate with a hole 
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We determine exactly the non-retarded Casimir-Polder interaction between a neutral but polar- 
izable particle and a perfectly reflecting sheet containing a circular hole. The calculation reveals a 
strong dependence of the interaction on the orientation of the particle's electric dipole moment with 
respect to the surface. For a dipole moment that is polarized perpendicularly to the surface the 
interaction potential has two saddle points lying above and below the plane of the surface, on a line 
through the centre of the aperture. For a dipole moment that is polarized parallel to the surface 
there is only one saddle point right in the middle of the aperture. Provided the particles' motion 
could be confined to a line through the middle of the aperture, this effect could potentially be used 
for population-sensitive trapping of particles. 



It is well-known that the interaction of an atom with a 
surface depends on the orientation of the atomic dipole 
moment with respect to the surface. For a ground-state 
atom and a flat, perfectly-reflecting surface this difference 
manifests itself only in marginally different pre-factors in 
the same kind of analytical expression [Ij. However, for 
more complicated structures such effects may be more 
pronounced. Recently we have shown that the large- 
distance behaviour of the interaction between an atom a 
cylindrical wire may display either logarithmic or power- 
like behaviour depending on how the atomic dipole is 
oriented with respect to the wire [2l|3]. Furthermore, it 
has also turned out that there is no Casimir-Polder at- 
traction for an atom placed exactly above the edge of 
a semi-infinite sheet of a conductor provided the atomic 
dipole moment is perpendicular to the material's surface 
[21 13] . The question arises as to whether this conclusion 
would continue to hold if we were to add to the system 
some other (infinitely) thin, flat conductors lying in the 
same plane. In particular, would the energy shift of a 
z-polarized dipole vanish when the dipole is located in a 
circular aperture of a perfectly reflecting sheet? (cf. Fig. 
[T]). This is interesting because we know that the atom- 
plate interaction energy at very large distances, where 
the presence of the aperture can be neglected, approaches 
zero from below [1 . Thus, the vanishing of the energy 
shift at z = within the aperture would imply that the 
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FIG. 1. (Color online) A dipole at a position ro in the vicinity 
of the perfectly reflecting sheet with a circular aperture of 
diameter d. 



potential calculated along any path that is perpendicu- 
lar to the plate and falls within the aperture features a 
minimum at some point [4^ . Such a minimum would im- 
ply some intricate behaviour of the Casimir-Polder force 
acting on a z-polarized dipole close to plate with a hole. 
Therefore, it seems worthwhile to investigate the problem 
in detail. To capture the essential features of the prob- 
lem we consider the non-retarded limit of the Casimir- 
Polder interaction, that is, we assume that the atom- 
surface separation is much smaller than the wavelength 
of the strongest dipole transition within the atom. This 
is the regime of the greatest practical interest, especially 
in the context of cold-atom experiments whose aim is to 
trap and guide atoms very close to surfaces. The non- 
retarded limit of the Casimir-Polder interaction between 
a neutral atom and a perfectly reflecting microstructure 
can be worked out by methods of electrostatics. Envis- 
aging the atom to be a point-like electric dipole located 
at ro one can show that the interaction energy with a 
nearby surface can be written in rectangular coordinates 
as % 
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where Gh(^^^') is the homogeneous part of the electro- 
static potential at r due to a point charge at r^ The sum 
runs over the three components of the electric dipole mo- 
ment operator /i^. The potential G(r,r') due to a unit 
point charge at satisfies 

- V'G(r,rO = (5(^Hr-rO (2) 

with appropriate boundary conditions. In the case of a 
perfect reflector, G(r, r') would have to vanish for any r 
on its surface. The homogeneous part of the potential 
is the difference between G(r,r') and the potential of 
the same point charge in free space, i.e. without any 
boundary conditions, 

G^(r,r') = G(r,r')-^-^. (3) 
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So far the problem is entirely classical and does not in- 
volve quantum electrodynamics. The quantum proper- 
ties of an atom are accounted for in the expectation 
values of the electric dipole moment operator = 
where \j) denotes the atomic state, not neces- 
sarily its ground state. The difficulty of determining the 
energy shift ([T]) lies in calculating the potential G(r, r') 
of a point charge, also called a Green's function, for the 
geometry of interest. In our case we need the solution of 
([2| for point charge placed close to a plate with a circular 
hole of diameter as shown in Fig. [l] In principle, this 
Green's function G(r, r') may be obtained by an expan- 
sion in terms of the eigenfunctions of the wave operator. 
For that the boundary conditions need to be imposed in 
oblate spheroidal coordinate system, where a plate with a 
hole is a surface of just one coordinate being constant, see 
e.g. [5], but then the eigenfunctions of the wave opera- 
tor are oblate spheroidal wave-functions whose handling, 
both analytically and numerically, is highly complicated 
[6l [7]. However, the required Green's function may be 
obtained by a far simpler method, namely, using a coor- 
dinate transformation discovered by Lord Kelvin [8 , the 
so-called Kelvin inversion. The Kelvin inversion is a non- 
linear coordinate transformation, a reflection of space in 
a sphere of radius S and centred at s, and is defined as 
follows: 



r\r] = 



r(r-s) + s. 



(4) 



The geometrical properties of this transformation are 
most clearly presented in [9]. The Green's function of the 
Poisson equation G(r, r') is effectively a function of two 
variables: the observation point r and the point where 
the source is located r^ To apply the transformation Q 
to the Green's function G(r, r') we apply it to both of its 
arguments. The transformation can be used to generate 
new solutions of the Poisson equation from the known 
ones. This is possible because we have [10] 



(5(^)(r-rO, (5) 



that is, the transformed Green's function G (T[r], T[r']) 
augmented by an appropriate pre-factor also satisfies the 
Poisson equation [11 . Thus, the transformation 
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r — s r' — s 



rG{T[r],T[r'])^G{r,r') (6) 



generates a new Green's function of the Poisson equation 
G(r, r') from an already known solution G(r, r'). The 
crucial point here is that, if the Green's function G(r, r') 
vanishes on some surface a, then the transformed Green's 
function G(r, r') vanishes on the transformed surface a, 
so that, by adjusting the position and the radius of the 
inversion sphere, one is able to obtain Green's functions 
for perfect reflectors of various shapes [TJl [13] . Here we 




FIG. 2. (Color online) Green's function ([7| corresponds to 
a half-plane geometry when the perfect reflector occupies the 
2/ = n X > plane (left). Application of the Kelvin inversion 
([6| with s = {d, 0, 0) and S = d yields the Green's function 
for a plate with a hole (right). 



follow this idea to generate the Green's function for a 
plate with a hole from the known case of the semi-infinite 
half-plane. The Green's function for a perfectly-reflecting 
half-plane expressed in cylindrical coordinates has been 
calculated in |2]: 
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The formula is valid for r and r' lying on the same side 
of the half-plane. The abbreviations in Eq. ^ are: 



2pp' cos(0 ^ (/)') + (z — 2: 



/\2 



= v'2ppqi + cos(0T^O]. e = sgn[sin(0 + 

The Green's function ^ satisfies equation ^ and van- 
ishes on the surface (j) = which corresponds to the 
surface y = D x > in Cartesian coordinates, 
as illustrated in Fig. [2] (left). In order to obtain the 
Green's function for a plate with a circular hole we ap- 
ply the transformation ([6| with the centre of inversion 
located at the point s = (d, 0, 0) and with its radius 
being S = d. For this particular choice of s and 
the transformation maps the half-plane <j = {r G : 
y = n x>0}toa sheet with a circular hole: 
a = {r eR^ :y = H (x - d/2)^ ^ > cf. 
Fig. [2] (right). Naturally we wish the coordinate system 
to be placed at the centre of the aperture in such a way 
that the z-axis is perpendicular to the sheet, as shown 
m Fig. [1] Thus we shift and rotate the axes according 
to X ^ d/2 -\- y, y ^ z ^ x. It is straightforward 
to check that the procedure described above yields the 
transformed Green's function in the form as Eq. ^ but 
with 
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and e = sgn[z(p'2 ^ _ ^2/4) ^ ^/(^2 ^ ^2 _ ^2/4)]^ ^he 
resulting Green's function is valid for r and lying on 
the same side of the sheet. It satisfies Poisson's equation 
([2| and vanishes on the surface of the conductor, i.e. for 
ah r with z = fl x'^ ^ y'^ > (P / 4. To see this it 
is instructive to go into spheroidal coordinates [5 where 
for ^, T], T]' all positive we get 



- - T ^0' + 2(1 + T 
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The knowledge of the Green's function now puts us into 
the position to calculate the atomic energy shift ([T]), 
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where we have defined 
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(12) 



The result (|8| together with Eqs. ([9|-([TT]) is, as far as 

electrostatics is concerned, exact. It applies to atoms 
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FIG. 3. (Color online) Interaction potential between a plate 
with a hole of diameter d and a dipole polarized in the p- 
direction, cf. Fig[l] The inset illustrates the plane on which 
Ep is evaluated and the orientation of the dipole. 
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FIG. 4. (Color online) Interaction potential between a plate 
with a hole of diameter d and a dipole polarized in the z- 
direction, cf. Fig^ The inset illustrates the plane on which 
is evaluated and the orientation of the dipole. 



whose distance from the surface is much smaller than the 
wavelength of the dominant dipole transition. The sim- 
plicity of Eqs. ([9|-([TT]) allows us to visualize the Casimir- 
Polder potential felt by an atom. As anticipated, for the 
z-polarized atomic dipole, the shift vanishes in the plane 
of the sheet, cf. Fig. [4j which can be confirmed also ana- 
lytically. The potential evaluated along lines of constant 
p < d/2 indeed has a minimum above and below the aper- 
ture. Interestingly, there is in principle a region where 
a z-polarized dipole is repelled from the microstructure, 
as has been concluded on the basis of numerical calcu- 
lations along the line p = [4 . One of the benefits of 
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FIG. 5. (Color online) The direction of the Casimir-Polder 
force acting on a z-polarized atom interacting with a plate 
with a hole of diameter d (top). The (red) dot shows the 
location of the saddle point at ^ ~ 0.371 Id. The direction 
of the Casimir-Polder force acting on isotropically polarized 
atom does not feature unusual patterns (bottom). 



our exact analytical result for the complete p-z plane is 
to show that motion along the line p = is unstable and 
upon the slightest perturbation the dipole is attracted 
towards the edge of the aperture, cf. Fig. |5j Thus, it 
is more appropriate to speak of attraction along unusual 
paths rather than of repulsion. The location of the sad- 
dle points and the magnitude of the shifts there can be 
estimated to be 



^ ±0.3711d, E,{p = 0,z= \zs\,d) 



4.0622 



(13) 



The potential experienced by the atom polarized in the p 
direction does not display the same peculiarities as that 
of a z-polarized one. We plot it in Fig. [3] and stress that 
the potential for the (/)-polarized atom looks very similar. 
In fact, both functions, Ep and S^, coincide at p = and 
feature a saddle point only at the centre of the aperture 
where we have 



^pAp = 0,z = 0,d) 



21.3333 



(14) 



Comparing Fig. [sjwith Fig. |4j together with Eqs. (13)- 
(14), we note that although the potential strongly de- 



pends on the polarization of the dipole the characteristic 
features of Sp^^ are by far more pronounced than those 
of Ez. Thus, for ground-state atoms with isotropic polar- 
izabilities the Casimir-Polder force will not display any 
unusual behaviour, cf. Fig. [Sj However, for cold polar 



molecules that by some other means are being confined 
to move along the z axis, one could in principle use this 
system for population-sensitive trapping: molecules with 
dipoles along z would trap above and below the plane 
of the aperture, whereas molecules with dipoles paral- 
lel to this plane would trap right at the centre of the 
aperture. The only drawback to bear in mind is that the 
relative depth of the minima of the two potentials is fixed 
to about 5.25 independent of the size of the aperture, so 
that the population of molecules trapped with dipoles 
along z would be smaller. 

It is worth noting that a similar polarization depen- 
dence of the dipole-surface interaction would be expected 
if one were to consider a sheet of a conductor with a 
long slit rather than a circular aperture. Such a geome- 
try might lend itself to polarization-sensitive guiding and 
scattering of cold polar molecules. 

The results presented here do not rely on the optical 
properties of the surfaces but rather on their geometry. 
Thus the perfect-reflector model considered here is fully 
appropriate. In practice the most important shortcoming 
of the model is the assumption of zero thickness of the 
conductor, which is the reason for the z-polarized dipole 
not to couple to sheet when it is located within the aper- 
ture. For cold atoms this condition is difficult to satisfy 
in practice, even in the era of graphene, and in reality we 
should expect the thickness of the material to come into 
play and alter the interaction. Numerical calculations 
have given estimates of how much the potential wells are 
being reduced by finite thickness of the material [4 . 
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